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This article is a pedagogical review of theoretical studies of noncentrosymmetric supercon- 
ductors with particular emphasis on the role played by electron correlation, which is important 
for heavy fermion systems. We survey unique properties of parity- violated superconductivity 
such as the admixture of spin singlet and triplet states, unusual paramagnetism, large Pauli lim- 
iting fields, magnetoelectric effects, the helical vortex phase, and the anomalous Hall effect. It 
is pointed out that these remarkable features are strengthened by a strong electron correlation 
effect, and thus are prominent in heavy fermion superconductors without inversion symmetry. 
We also discuss possible pairing states realized in the heavy fermion system CePtsSi. 
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1. Introduction 

In all superconductors, the symmetries of systems such 
as time-reversal symmetry and inversion symmetry im- 
pose important constraint on pairing states. For instance, 
when a system has an inversion center, Cooper pairs are 
classified into a spin-singlet state or a spin-triplet state 
according to parity. Conversely, in the case without inver- 
sion symmetry, in which an asymmetric potential gradi- 
ent yields a spin-orbit (SO) interaction that breaks par- 
ity, the admixture of a spin singlet state and a spin triplet 
state is possible, and also, various exotic electromagnetic 
properties may result. Such parity- violated superconduc- 
tivity is the main topic of this article. Pioneering theoret- 
ical studies on this subject were initiated by Bulaevskii ct 
al. decades ago, and developed by Edelstein in 1989. 1-4 
Edelstein pointed out that, in addition to the singlet- 
triplet mixing of pairing states, the parity-breaking SO 
interaction also gives rise to unusual electromagnetism 
such as van-Vleck-like spin susceptibility analogous to 
van-Vleck orbital susceptibility, and magnetoelectric ef- 
fects, i.e., supercurrent induced by the Zccman field, and 
conversely, magnetization induced by supercurrent flow. 
The latter effects are due to the nontrivial coupling be- 
tween the spin and charge degrees of freedom raised by 
a parity-breaking SO interaction. These intriguing phe- 
nomena were also investigated in more detail by Gor'kov 
and Rashba, and Yip. 5 ' 6 However, experimental tests on 
these interesting theoretical predictions in real materi- 
als have not been performed for decades since the first 
prediction made by Edelstein. 1 

In 2004, Bauer et al. discovered a bulk superconduc- 
tor without inversion symmetry, CePtaSi, 7 and stimu- 
lated more extensive experimental and theoretical stud- 
ies on this subject. 8-27 In particular, the ETH group and 
Agterberg revealed some novel aspects of electromag- 
netic properties such as large Pauli limiting fields, and 
a possible realization of the helical vortex state in mag- 
netic fields. 15 ' 19 The helical vortex state is a kind of the 
Fulde-Ferrel state in which Cooper pairs have a center of 
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mass momentum, and the phase of the order parameter 
is spatially modulated. This Fulde-Ferrel-like state was 
also independently predicted by Samokhin on the basis 
of the general GL theory with strong SO coupling. 16 

After this breakthrough, other novel heavy fermion su- 
perconductors without inversion symmetry such as Ulr, 
CeRhSi3, and CelrSia have been discovered. 28-30 Also, 
in non- /-electron systems, new noncentrosymmetric su- 
perconductors such as Cd2Re207, Li 2 Pd3B, and Li 2 PtaB 
have been discovered. 3 1-34 An intriguing feature of these 
systems is that according to LDA calculations and de 
Haas-van Alphen experiments, the SO splitting of the 
Fermi surfaces is much larger than the superconducting 
gap. 17, 35-37 For instance, in the case of CePtaSi, the mag- 
nitude of the SO splitting is nearly equal to 10 percent of 
the Fermi energy. 17 ' 35 This large SO splitting indicates 
the important role of parity violation in these supercon- 
ductors. 

In this article, we survey theoretical studies of noncen- 
trosymmetric superconductors with particular emphasis 
on the role played by electron correlation effects and the 
application to heavy fermion systems. In heavy fermion 
systems where /-electrons arc itinerant, there exists a 
strong spherically-symmetric SO interaction in addition 
to the parity-breaking SO interaction. Thus, we need to 
clarify the relation between these different types of SO 
interaction. The parity-breaking SO interaction due to an 
asymmetric potential gradient W is generally expressed 
as (fe x VV) • cr, where fc is the momentum, and <x is the 
Pauli matrix. This SO interaction is expressed in terms 
of electron spins. In the case of CePtaSi, because of the 
strong symmetric SO interaction, the ground state of the 
/-electron level is the Tr Kramers doublet, the basis of 
which is, 9 



0/61 -3/2), 



5/2) 



(1) 



Transforming the basis of spin to the basis \a), \a), we 
find that the parity-breaking SO interaction is rewrit- 
ten as (5/12) (fc x VV) • c, where the Pauli matrix <r 
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is denned in the pseudospin space spanned by \a), \a). 
That is, the matrix structure is not affected, and onfy 
the constant prefactor is changed. Thus, even in the case 
of heavy fermion systems, as long as the Bloch states of /- 
electrons are labeled by pseudospins that constitute the 
Kramers doublet, we can use the same expression of the 
parity-breaking SO interaction as that for electron spins. 
Throughout this paper, when we consider heavy fermion 
systems, we use the terms "spin singlet state" and "spin 
triplet state" to indicate the pseudospin singlet state and 
the pseudospin triplet state, respectively. In the case that 
upper levels of /-electrons are not negligible, we need a 
more elaborate treatment of the parity-breaking SO in- 
teraction. 

The organization of this paper is as follows: In §2, we 
present the microscopic description of the superconduct- 
ing state without inversion symmetry which takes elec- 
tron correlation effects into account. This section is a bit 
technical. Readers who are not interested in the technical 
details may proceed to the next sections. In §3, we argue 
the pairing state on the basis of the BCS weak-coupling 
approximation, focusing on the admixture of spin singlet 
and triplet states. In §4, we study unique electromag- 
netic properties of noncentrosymmetric superconductors. 
In §5, we provide a discussion on possible pairing states 
realized in CePtsSi. In §6, we give concluding remarks. 

2. Fermi Liquid Theory for Superconductors 
without Inversion Symmetry 

In this section, we present the basic formulation for 
the BCS theory of noncentrosymmetric superconductors, 
which includes electron correlation effects as the Fermi 
liquid corrections. 22 ' 27 We are mainly concerned with 
electromagnetic properties. As is shown below, unique 
features of magnetism in inversion-symmetry-broken sys- 
tems appear as a result of paramagnetic effects rather 
than diamagnetic effects. Thus we take into account 
the Zeeman coupling with an external magnetic field 
H = (H x , H y , H z ), neglecting the orbital diamagnetic 
effect for a while. Thus, our analysis is based on the fol- 
lowing simple Hamiltonian, 

(2) 
(3) 

(4) 
(5) 



n 

Ho 



Hso 



— Ho + Hso + ^Zeeman, 

= ^ ekcjck + U~y] nana 



H 



Zeeman 



k 

c\n B a ■ Hc k . 



k 



where c\. = (c| fe , c{ fe ) is the two-component spinor field 
for an electron with the spin f, J., and the momentum 
k. n a i — c^Cai is the number density operator at the 
site i. The components of <x = (a x , <7 y ,<j z ) are the Pauli 
matrix. The vector Co(k) = (£ox> £-0y> A)z) m the spin- 
orbit interaction term Hso obeys Co(—k) = —£ (k) 7 
breaking inversion symmetry. The explicit form of Co(k) 
is determined by details of the crystal structure. In the 
case of the tetragonal structure with an asymmetric po- 
tential gradient along the z-axis, Co{k) = (k y ,—k x ,0) 



for a small fc, which is the Rashba interaction. 38 In the 
case of cubic structures with the point group symmetry 
T, such as the zinc-blende, Co{k) = (k x (k y — kl) 7 k y (k z — 
k x ),k z (k x — ky)), which is related to the Dresselhaus 
interaction. 39 For simplicity, we assume that the g- 
f actor is equal to 2, and spin s = 1/2. As a matter of 
fact, for heavy fermion systems the g- factor may gener- 
ally be different from 2 because of the strong symmetric 
SO interaction. 

Electron correlation effects, which are important for 
heavy fermion systems, are incorporated by on-site 
Coulomb repulsion U. We also assume that, in addition 
to the on-site repulsion, there is an effective pairing in- 
teraction with an angular momentum I > 1, which may 
stem from the on-site Coulomb interaction in (2) or from 
any other factors not included in the Hamiltonian (2). 

To calculate physical quantities microscopically, we 
utilize the single-particle Green function, the inverse of 
which is, in the conventional Nambu representation, de- 
fined as 



9 {P) ~\ -At(p) 



-A(p) 



(6) 



where p = (is n , k) and 

H(p) = e k - n + aC (k) ■ <x + £{p) + /i B er • H, (7) 

with fi being the chemical potential. The self-energy ma- 
trix £ consists of both diagonal and off-diagonal compo- 
nents, and expressed as 

y - ( E Tt(rt s Tl(rt 
S iT (p) E u (p) 

= S + S • cr, (8) 

where S = (S x , £„, E z ) with S = (S TT + E u )/2, Y> x = 
(E iT + E n )/2, Y, y = (E iT - £ n )/2i, and E z = (E TT - 
S U )/2- 

The SO interaction aCo(k) ■ cr splits the Fermi sur- 
face into two parts as shown in Fig. 1. The SO split- 
ted two bands are derived by diagonalizing ie n — H(p) 
and ie n + H f (—p) in Q~ 1 (p). The diagonalization is a bit 
tricky because S(p) is generally non-Hcrmitian. In the 
absence of the magnetic field, we can verify that E x (p) 
and E y (p) are real quantities, and thus is Hermi- 

tian. 27 Also, as long as the linear response against the 
magnetic field is concerned, the imaginary part of H x (p) 
and Ej/(p) are negligible compared to the real parts even 
for H 7^ 0. Then, the diagonalization is achieved by using 
the unitary transformation A(p)G~ l '(p).A' (p) with 



A(p) 



U(p) 








dip) 



£±(p) 



v±(p) 



l 

71 



i± 



4(p) 



C x (p)±iC y {p) 
yf 'C x {pf + C y {pf 



(9) 
(10) 

(11) 
(12) 
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Here, 

C{p) = (C x {p),C v {p),C x {p)) = C {k) 



Mb 



H 



and \\C(p)\\ = 

In the absence of the magnetic field H = 0, the su- 
perconducting gap function A(p) in eq. (6) is also diag- 
onalizcd by A if and only if the gap function has the 
form, 1 ' 15 ' 22 

A(p) = A s (k)ia y + A t (k)C(p) ■ tri(j y 

That is, when the d-vector of the triplet gap satisfies 
d oc C(p), there is no Cooper pair between electrons 
on the different Fermi surfaces splitted by the SO in- 
teraction. Conversely, unless d oc C(p), pairing corre- 
lation between the SO splitted two bands is induced, 
which yields pair-breaking effects, (see Fig.l.) Accord- 
ing to weak-coupling calculations performed by Frigeri 
et al., 15 the d- vector parallel to C(p) gives the highest 
T c , as long as the pairing interaction stabilizes the gap 
function with the same momentum dependence as that of 
C(p). In general, a finite magnetic field gives rise to pair- 
ing correlations between the different Fermi surfaces, and 
thus the gap A(p) can not be diagonalized by A. How- 
ever, within the linear response theory, which is adequate 
for the following argument, the correlation between the 
different Fermi surfaces is negligible, and the inverse of 
eq.(6) is given by, 



Q{p) 



G(p) F(p) 
F\p) -G\-p) 



where 



T=±l 



Hp) 



El + rt(p) ■ er . 
ta y F T (p), 



r=±l 



and 



Grip) = 
Fr( P ) 



z kT {ie + e* kT ) 



[ie + ij kT sgne) 2 - E\ T ' 

ZkrAkr 



{ie + i7 fer sgne) 2 - E\ T ' 
Here, E kr = ^e* 2 + A 2 r , A kr = z kr A kr with 

A fcT = A s (fc)+r|£(fc)|A t (fc), (20) 
and e* kT is the quasiparticle energy determined by using 
e* kT + m - e k - S «, k) - ra\\C(e* kT , k)\\ = 0. (21) 
7fe T is the quasi-particle damping given by 



In the derivation of the above equations, we assumed 
that the quasiparticle approximation is applicable. As 
-S(p). (1^11 be seen in the following, when the SO splitting is 
much larger than the superconducting gap, some elec- 
tromagnetic properties are dominated by electrons far 
from the Fermi surfaces. In such cases, the quasiparticle 
approximation is not generally applicable, and one needs 
to take into account the full-energy dependence of the 
self-energy. 

(14) 3. Pairing State — Admixture of Spin Singlet 
State and Spin Triplet State 

One of the most remarkable features of noncentrosym- 
metric superconductors is the admixture of spin singlet 
and spin triplet states. 1 ' 5,15 This phenomenon is concep- 
tually understood as follows. When a system has a crystal 
structure without an inversion center, an asymmetric po- 
tential gradient VT^ exists, causing the spin-orbit (SO) 
interaction (k x W) ■ cr, which breaks parity as well 
as spin inversion symmetry. The asymmetric SO interac- 
tion splits the Fermi surfaces into two parts. We show a 
simple example in Fig.l in the case of W \\ n = (001), 
which corresponds to the Rashba-type SO interaction. In 
Fig. 1, the spin quantization axis is parallel to kpx W. 
Then, on one of the SO splitted Fermi surfaces, a Cooper 
pair consists of an electron with the momentum k and 
spin — », and an electron with the momentum — k and 
spin <— ; i.e. \k — »)| — k «— ). Since the counterpart of this 
state, | <— ) | — >), is formed on the other Fermi surface, the 

(15) superposition of these states is impossible. The pairing 
state | — >) | «— ) is an admixture of a spin singlet state and 
a spin triplet state with an in-plane spin projection equal 
to zero, i.e., | -)| <-> = [(| -)| <-> - | <->| -)) + (| - 

( 16 ) ) | <— ) + 1 <-)| — >))]/2. This means that for the spin quan- 
tization axis parallel to the z-axis, the spin triplet com- 
ponent has the nonzero spin projection, S z = ±1. (see 

(17) Fig. 2.) More precisely, the spin part of the triplet pairs is 
determined by the parity-breaking SO interaction, ensur- 
ing that the pairing correlation between different Fermi 
surfaces, which yields pair breaking effects, is suppressed, 

(18) as pointed out in the previous section. 

To obtain further insight into the pairing state, we 
switch to a microscopic argument based on the BCS the- 

(19) ory. The superconducting gap function A Q( 3 and the tran- 
sition temperature are determined by the self-consistent 
gap equation, 



A a p = Tj2nr af3 (p, P ')F(p% 



(24) 



7fer = z kT (lm^(e, k) + aRet ■ lmS R (e, k))\ e 



(22) 



where ReC — ReC R (e, k)/\ReC R (e, k)\ and quantities 
with the superscript R indicate retarded functions ob- 
tained by the analytic continuation ie n — > e + iS (5 > 0). 
The mass renormalization factor z kr is defined as 

-l 



Z kT 



1 dXo(p) ra d\\C(p)\\ 



d(is n ) 



d(ie n ) 



(23) 



where F(p) is the anomalous Green function defined by 
cq. (17). We have introduced the four-point vertex func- 
tion matrix {T a P (p,p')} 7 s- This pairing interaction con- 
sists of the parity-conserving part r con and the parity- 
nonconserving part f noncon , i.e., f = f con + r noncon . The 
former can be decomposed into spin singlet and triplet 
channels. Let us assume the following form of the pairing 
interaction. 

{T^ n ( P ,p')hs = rgLM)i(<r v )af>i((rvh6 
+rW,(p,j/)£(p) • (<ria v ) af3 C(p') ■ (aia v ) l5 (25) 
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Fig. 1. Two Fermi surfaces splitted by the Rashba SO interaction. 
On each Fermi surface, spins are parallel to fc X VV, and have 
a definite eigen value of (fc X n) ■ <r/|fc|; i.e. +1 or -f. When 
d oc jCo(k), Cooper pairs are formed on each Fermi surfaces 
(black and gray dotted thin arrows). When the d- vector is not 
parallel to Co(k), the pairing correlation between the different 
Fermi surfaces exists (black dotted thick arrows), which gives 
rise to pair-breaking effects. 




H 

or 

III 



Spin Triplet 

Fig. 2. Schematic picture of parity-violated Cooper pairs. The 
left-hand side represents the Cooper pair | — >)| <— ) with the 
spin quantization axis parallel to fc X n. The right-hand side 
represents the sum of a spin singlet state and a spin triplet state 
with S z = ±1 for the spin quantization axis parallel to the 2-axis. 



The form of the triplet part is chosen to prevent the pair- 
ing between the different Fermi surfaces due to the mis- 
match between the d- vector and the momentum depen- 
dence of the SO interaction, as pointed out in the previ- 
ous section. 1 ' 15 ' 22 The four-point vertices Tcon {v = s, t) 
are expanded in terms of the basis of the irreducible rep- 
resentations of the point group of the system. Generally, 
the parity-nonconserving part r noncon is nonzero, and ex- 
pressed as 



expressed by the same basis function of the irreducible 
representation. The superconducting state realized is the 
mixture of the spin singlet and triplet states. 1,5 In this 
case, the triplet pairing state should be in a higher angu- 
lar momentum state than the singlet pairing state, and 
the possible pairing state is s + p or d + / or g + h, and 
so forth. 

To examine to what extent the admixture is realized, 
we carry out a simple model calculation of the transition 
temperature T c and the gap function A. We consider a 
two-dimensional model with an isotropic Fermi surface 
and the Rashba interaction Co(k) — (k y ,—k x ,0). We 
also assume the existence of a pairing interaction in the 
s-wave and the p-wave channels expressed as 



s c lt c -fel c - fe 'i Cfc 'T 



+Wp[(ifc + )(ifc+)4 T cl fcT c_ fc ' T Cfc' T 



+ {ik-)(ik'_ 



)cl l cl kl c- kll c kn ], (27) 



with k± = (k x 
tion processes, 
tion is 



± iky)/\k\, and ignore any other interac- 
Then, the weak-coupling BCS gap equa- 



A T = ^2 9tt> ^2 A T . 



T'=± 



tanh( ^y' 

2-Efcr' 



(r = ±) 



(28) 



where A T = A s + rA t with A s and A t the gap functions 
for the spin singlet and spin triplet states, respectively. 
9tt> = ~(u s + TT f u t )/2 and E kr = ^e z kT + A 2 . The 
transition temperature T c is easily calculated by using 
the linearized gap equation. In Fig. 3 (a), we show the re- 
sults of T c plotted as a function of u s for the attractive p- 
wave interaction with u p — -0.15.Ef, and Ego/ Ep = 0.1 
or Eso/Ep = 0.05 with Eso the magnitude of the SO 
splitting. For the attractive s-wave interaction u s < 0, 
T c increases, as |u s | changes from to indicating 
that the s-wave pairing mixes with the p-wave pairing 
cooperatively. It is noted that for u s w u p , T c is substan- 
tially increased by the admixture. On the other hand, for 
the repulsive s-wave interaction u s > 0, T c is decreased 
by the admixture. The decrease in T c is at most of or- 
der Eso I Ef for all positive values of u s including an 
extreme case of u s /\u p \ oo. To observe the admixture 
more directly, we show in Fig. 3(b) the calculated results 
of the gap function for the spin triplet pair A t and the 
spin singlet pair A s in the case of Eso/Ep = 0.1 at zero 
temperature. For the attractive u s < 0, the mixed s-wave 
gap grows significantly, as |u s | increases. The maximum 
ratio of A s /A t for u p < u s < is nearly equal to 0.32. In 
contrast to the substantial admixture for u s < 0, A s /A p 
is at most of order Eso/Ep = 0.1 in the case of the 
repulsive u s > 0. We found that even in the limit of 



» i . . i f u s — > oo, there exists a solution of the gap equation: 

{ r noncon(>>P')}7<5 = t(s (p,p')i{°v)<*p£{p') ' (<™<7 y ) 7 a Af = q.002498, A s = 0.0001468. Thus, even when there 

+T {ts \p,p')C{p) ■ (cTia y ) a pi{cT y ) 1 8. (26) is a strong on-site Coulomb repulsion, as in the case of 

heavy fermion systems, the s+p wave state may be pos- 
and r(* s ) are also classified in accordance with the s ible, as long as Eso/Ep is sufficietly smaller than unity. 

In the above argument, it is postulated that the d- 
vector of the triplet pairing can be compatible with the 
momentum dependence of the SO interaction Co(k) to 
prevent the pairing between the different Fermi surfaces. 



p(ai 

irreducible representations. The highest T c is achieved 
by r( v '(p,p') (y — s,t, st,ts) which belongs to a unique 
irreducible representation, provided that there is no acci- 
dental degeneracy. Thus, A s (p) and A t (p) in eq.(14) are 
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Fig. 3. (a) T c vs u s . (b) A s and At vs u s at zero temperature 
for E S o/E F = 0.1. 



However, this condition is not generally guaranteed. It 
may be possible that for some specific electronic struc- 
tures, the strongest attractive interaction in the triplet 
channel has a momentum dependence inconsistent with 
Co(k), and thus the gap function A(k) cannot be diago- 
nalized by the unitary transformation Z7(fe)A(fc)t/*(— k). 
The realization of this pairing state depends on the com- 
petition between the strength of the attractive interac- 
tion and the pair-breaking effect caused by the pairing 
between the different Fermi surfaces. 

4. Electromagnetic Properties 

4-1 Paramagnetism 

The parity-breaking SO interaction yields unique para- 
magnetic properties. For example, the transition between 
the two SO splitted Fermi surfaces gives rise to a van- 
Vleck-like spin susceptibility in addition to the Pauli spin 
susceptibility 1,5 ' 15 ' 22 The van-Vleck-like term is analo- 
gous to the conventional van-Vleck orbital susceptibil- 
ity. However, there are some crucial differences between 
them, as argued below. An important feature of the van- 
Vleck-like term is that when the SO splitting Eso = a\C\ 
is much larger than the superconducting gap A, as in the 
case of CePtaSi, CeRhSi3, and CeIrSi3, this term is not 
strongly affected by superconducting transition, 1 ' 5,15 in 
contrast to the Pauli term which decreases below T c for 
spin singlet pairing states. 

Here, we first consider the Rashba SO interaction rele- 
vant to tetragonal systems, Co(k) — (t 0y , —t 0x , 0), where 
tov iy — x,y) transforms like k v . For small |fe|, t^^ = k^. 
In the case of CePtaSi, CeRhSi3, and CeIrSi3, the crystal 
structures of which have d v symmetry, the z-component 
of jCo is generally nonzero, but given by k x k y k z (k 2 — /c 2 ) 
for small \k\, as elucidated by Samokhin. 21 However, to 
highlight the effects of the van-Vleck-like term, we ex- 
ploit the Rashba type interaction for simplicity. From 



(a) 



(b) 



H, 





Fig. 4. Effects of the Zeeman magnetic field on the changes of the 
Fermi surfaces, (a) For magnetic field parallel to z-axis. Arrows 
on the Fermi surfaces indicate spins. The length of the arrows 
represents the magnitude of magnetic moments, (b) For in-plane 
magnetic field. 



cq. (21), the single-particle energy in the absence of the 
Coulomb interaction for H = (H x , 0, H z ) is given by 

e k ± =e k ±ay/ (t 0y - fi B H x ) 2 + t\ x + ^i B H 2 z . (29) 

It is seen from eq. (29) that magnetic fields parallel to the 
z-axis merely change the magnitude of the SO splitting 
of the Fermi surfaces, while magnetic fields parallel to 
the rry-plane deform the Fermi surfaces into asymmetric 
shapes, as depicted in Fig. 4. These effects indicate that 
magnetic fields parallel to the z-axis is dominated by the 
van-Vleck-like susceptibility, while the in-plane magnetic 
fields induce responses from electron spins corresponding 
to both the Pauli term and the van-Vleck-like term. 

According to the Fermi liquid theory presented in sec- 
tion 2, the exact expression of the spin susceptibility in 
the superconducting state for the Rashba interaction is 
obtained as 27 



,_. 2 \ ~> \ - tanh( 
Xzz{T) = A*b 2^ TZk ' 

k t=± 



2T 



EkT 



E kr +4+4- 



A fe+ A fe _ 



El + ~El_ 

for a magnetic field parallel to the z-axis, and 

Zk 



A sz (E kT ,k), (30) 
and 

t y A s T x (E kT ,k) 



T=± k 



+Mi£E^ tanh(fp) 



k t=± 



,*t + 4 + 4- 



E kT 
A fe+ A fc _ 



El 



■t x A s + x _(Ekr,k), (31) 



for an in-plane magnetic field. Here E kT and A kT are 
defined by the equations just below eq. (19) The three- 
point vertices A sz , A^_ , and K S T X are given by 



A"(p) = 1 - 



1 dX (p) 
Mb dH z 



1 dYl x ^ t x dY]y 



5S 



yU B dH x /j, b dH x fi B dH x 



^ B dH x fi B dH x ' 



(32) 



(33) 



(34) 
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Here, (t x (p),t„(p)) - i(p) = t(p)/\t(p)\, t(p) = n x £(p) 
with n = (001). Xzz(T) and the second term of Xxx(T) 
are the van-Vleck-like susceptibilities. Both the Pauli and 
van-Vleck-like terms are enhanced by electron correlation 
effects, because of the three-point vertices. 

It is noted that the van-Vleck-like term is quite dif- 
ferent from the usual van-Vleck orbital susceptibility, 
which is caused by symmetric SO interaction, for the fol- 
lowing reason. In contrast to the usual van-Vleck term, 
which is almost independent of temperatures, the van- 
Vleck-like term due to the parity-breaking SO interac- 
tion indeed depends on temperatures, similarly to the 
Pauli susceptibility, since the magnitude of the parity- 
breaking SO splitting is much smaller than the Fermi 
energy for typical heavy fermion systems. According to 
a simple model calculation based on a two-dimensional 
tight-binding model on a square lattice carried out by the 
author, both the Pauli and van-Vleck-like terms exhibit 
strong temperature dependence, as a result of the energy 
dependence of the density of states. 27 This property im- 
plies that it is almost impossible to distinguish between 
these two contributions experimentally by analyzing the 
temperature dependence of spin susceptibilities. Thus, 
we need to be careful about the interpretation of experi- 
mental data of spin susceptibilities in the superconduct- 
ing state, since they always include the van-Vleck-like 
term which is not strongly affected by superconducting 
transition when A <C -Esso- 
in the case of the Rashba SO interaction, the temper- 
ature dependence of spin susceptibilities below T c was 
calculated by Gorkov and Rashba, and Frigeri et al. for 
a model with a spherical Fermi surface. 5 ' 15 According 
to their results, for A <C Eso, Xxx at zero temperature 
is one half of the magnitude in the normal state, i.e., 
Xxx{Q) = Xxx(T c )/2, owing to the existence of the van- 
Vleck-like term. 

However, experimental observations seem not to be 
consistent with the theoretical prediction. According to 
the Knight shift measurements performed by Yogi et 
al. ' for CePt3Si, both Xxx and Xzz show no sig- 
nificant change even below T c . To explain the origin 
of the discrepancy, we note that the magnitude of the 
ratio Xxx(T = 0)/xxx(T = T c ) in the superconduct- 
ing state crucially depends on the details of the elec- 
tronic structure and electron correlation effects. For in- 
stance, when the density of states has a strong energy 
dependence, Xxx(T = 0)/xxx(T = T c ) may deviate 
from 1/2. To demonstrate this, we consider a Hubbard- 
type model on a square lattice with a kinetic energy 
Sk = —2{cosk x + cos k y ), and the Rashba term with 
C (k) = (sin sin ^,0). For this model, at half- 
filling, we calculate Xxx{0)/Xxx(T c ) in the case of A <C 
a\C\, handling electron correlation effects by perturba- 
tive expansion in terms of U up to the second order 
and using eq. (31). The results are shown in Fig. 5 As 
U/W increases, the van-Vleck-like susceptibility is more 
strongly enhanced by electron correlation than the Pauli 
term, and the ratio Xxx(T = 0)/xxx{T = T c ) approaches 
1.0. Although the strength of the electron-electron inter- 
action for which the ratio is nearly equal to 1.0 is consid- 
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Fig. 5. Xxx(T = 0)/xxx{T = T c ) in the superconducting state 
versus U/W with a = 1.0 calculated by second order perturba- 
tive expansion with respect to U for the 2D Hubbard-likc model. 



erably large, the results may be qualitatively valid, and 
can be improved by including higher order corrections, 
as long as the system is in the Fermi liquid state. As a 
matter of fact, CePt 3 Si is in the Fermi liquid state with 
a mass enhancement factor of order ~ 60, in which the 
spin susceptibility is also substantially enhanced. Thus, 
it is possible that, for CePtaSi, because of strong electron 
correlation and a particular electronic structure, Xxx{T) 
is dominated by the van-Vleck-like term, and thus is al- 
most independent of temperature below T c , consistent 
with the NMR measurements. 

In the case with cubic symmetry, the uniform spin sus- 
ceptibilities in the superconducting state are obtained in 
a similar manner. 27 The results are 

ZkT 



Xzz(T) = MbEE 



4Tcosh 2 (%f) 



Ap ub 



(E kT ,k) 



+MiEE^ tanh( ^ ) 



k r=± 



El 



E 2 k+ ~ El 



EkT 

Ak+Ak - A^( EkT ,k), 



(35) 



where 

A P ub (p) 



~ 2 i aso(p) 

£ z£xdY; x (p) C 0z t y dY ly (p) 



Mb 



dH z 



Mb 



dH z 



,(36) 



A<? h (p) = (£* +£*){! - 



1 dSuCp) 
dH z 



Mb 



+ 



C 0z £ x dY, x (p) | C 0z C y dY, y (p) 



MB 



dH z 



Mb dH; 

Because of the cubic symmetry, Xxx — Xyy — Xzz, and 
Xzz consists of both the Pauli term (the first term on the 
right-hand side of eq. (35)) and van-Vleck-like term (the 
second term of (35)). For systems with a spherical Fermi 
surface, the Pauli term is 1/3 of the total susceptibility. 
In this case, at zero temperature in the superconducting 
state, we have Xzz(0)/Xzz(T c ) = 2/3 for A < E so . 

4-2 Pauli limiting field 

In this section, we consider the Pauli depairing effect 
due to the Zeeman magnetic field. This effect in the 



J. Phys. Soc. Jpn. 



Full Paper 



7 



case without inversion symmetry was precisely studied 
by Frigeri et al. 15 The changes of the Fermi surfaces due 
to the Zeeman magnetic fields depicted in Fig. 4 imply 
that the Pauli depairing effect on Cooper pairs is also 
unique in noncentrosymmetric superconductors. 

We first consider the Rashba case with an asymmetric 
potential gradient along the z-axis. For magnetic fields 
parallel to the z-axis, Cooper pairs between electrons 
with the momenta k and — k are always possible, as 
shown in Fig. 4(a), implying the suppression of the Pauli 
depairing effect. More precisely, the Pauli limiting field 
depends on the spin state of Cooper pairs. For spin- 
triplet pairing states with the <i-vector determined by 
the Rashba interaction, d cx (to y , — tox, 0), the pairing 
state with S z = ±1 for the spin quantization axis par- 
allel to the z-axis is realized. Magnetic fields parallel to 
the z-axis do not give rise to the Pauli depairing effect 
on this state. On the other hand, for spin singlet states, 
the Pauli limiting field exists. However, it is strongly en- 
hanced by the parity-breaking SO interaction. To argue 
the Pauli depairing effect on the singlet state, we consider 
a simple case without electron correlation, and assume 
the attractive interaction —VwkW-k for a singlet chan- 
nel, where Wk is a structure factor for the orbital degrees 
of freedom. We also assume that Eso/Ep is sufficiently 
small, and that the admixture of triplet states is negligi- 
ble. Then, the gap equation for H z ^ is given by 



f 

V 



= EE 

k r=± 

+ EE 



a l*o| w l tami 



2T 



w\ii\Hl (kr tanh fff 



where 

Ckr = 1 + T 



a 2 \t a \ 2 + M |ff z 2 

v^wTj^m+T^HWk 



(38) 



(39) 



E kT =±[e\ + a 2 \t \ 2 + fi 2 B H 2 + A 2 

+2r ^s 2 +A 2 (e + -e) 2 } km 



and |io| : 



t()x + toy ^-k 



A Wk- It is evident that for 



HbH z <C ct|to | , the first term on the right-hand side of 
eq. (38) dominates over the second term, and thus the 
Pauli depairing effect is negligible. For general values of 
H z , the transition temperature T c is determined by using 
the linearized gap equation, ft is noted that the first term 
on the right-hand side of (38) in the linearized approxi- 
mation yields a log T singularity at low temperatures for 
any finite values of H z which satisfies [i^,H z < Ep, while 
the second term exhibits no logarithmic singularity. As a 
result, the finite transition temperature T c always exists 
for any finite magnetic field H z < Ep. This implies that 
the Pauli limiting field is infinitely large for sufficiently 
low temperatures. According to precise numerical cal- 
culations performed by Frigeri et al., the Pauli limiting 
field determined by the linearized gap equation obeys a 
concave curve on the H-T plane at sufficiently low tem- 
peratures. 15 ft is noteworthy that the recent experiments 
on CeRhSi3 support the existence of the large upper crit- 
ical fields at low temperatures exceeding a conventional 



Pauli limit, and exhibiting a concave behavior as a func- 
tion of T. 42 In these experiments, although the concave 
behavior is observed even near T c , it is expected that the 
inclusion of the strong-coupling effect enlarges the tem- 
perature region where the concave behavior appears up 
to higher temperatures close to T c . 

Also, the above results are easily obtained from a sim- 
ple thermodynamic argument. Let Xzz,n and Xz Z ,s be 
the spin susceptibilities in the normal and superconduct- 
ing states, respectively, for a magnetic field parallel to the 
z-axis. Then, the Pauli limiting field Hp at zero temper- 
ature is determined by equating the ground state energy 
in the normal state to that in the superconducting state 



Xzz,N 



Hi 



Xzz,sH 2 P N(Q)A 2 



(41) 



which leads to H P = A o ^/N(0)/(xzz,n - Xzz,s), where 
N(0) is the average density of states of the SO splitted 
two bands. We see from eq.(30) that when A <C Eso, 
Xzz,n~Xzz,s ~ Xzz,n(A /E so ) 2 < Xzz.n holds, leading 
to a large H P enhanced by the factor E$o/ A . More pre- 
cisely, the free energy difference between the supercon- 
ducting state and the normal state at any finite tempera- 
tures is calculated from O s -0„ = (dV' 1 /dA)A 2 dA. 
By carrying out numerical calculations, one can confirm 
that when the gap equation (38) gives a nonzero value 
of A , fl s — Q n < holds at least for A < -Eso, an d 
thus the superconducting state is the absolute minimum 
of the free energy; i.e. Hp determined by using the lin- 
earized gap equation coincides with that obtained from 
the thermodynamic consideration. 

For magnetic fields parallel to the xy-plane, the Fermi 
surfaces are deformed into asymmetric shapes as shown 
in Fig. 4(b), which gives rise to the Pauli depairing ef- 
fect for both spin singlet and spin triplet states. The 
Pauli depairing effect stems from the Pauli term of 
the spin susceptibility Xxx, i-e., the first term of eq. 
(31). When A < E so , and x™ n / lock - X™ n / lock « 

..van Vlcck / A / 771 \2 *y ...van Vlcck ...Pauli fV>i-. 

Xxx,n (Aq/Eso) < X xx ,n 1 Xxx.n nolds, the 
Pauli limiting field is H P w A N (0) / (xlxS ~ X™"#)- 
In the case with a spherical Fermi surface but with- 
out electron correlation, XxS"/v ^ s one na ^ °f tnc to " 
tal spin susceptibility. Thus, at zero temperature where 
Xxx^S = 0, Hp is enhanced by a factor of y/2, com- 
pared to the case with inversion symmetry in which the 
Pauli limiting field is ff^ v s y m - = A /V^- However, the 
real situation is more complicated. According to a more 
precise analysis done by Kaur et al. and Samokhin, the 
asymmetric deformation of the Fermi surface due to in- 
plane magnetic fields may stabilize Cooper pairs with a 
center of mass momentum, leading to the Fulde-Ferrel 
state. 16,19 We would like to discuss this point in section 
4.6. 

In the case with cubic symmetry, the spin suscep- 
tibility consists of both the Pauli and van Vleck-like 
terms, and an argument similar to that for the Rashba 
model with in-plane magnetic fields is applicable. When 
A <C Eso holds, the Pauli limiting field is given by 

Hp « Ao v /iV(0)/(xL aU A> i -xL au i i )- Here, X ^f {N) is the 
first term on the right-hand side of eq. (35) in the su- 
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perconducting (normal) state. In the case with a spher- 
ical Fermi surface without electron correlation, at zero 
temperature, Xzz at 1S °f the total susceptibility, as 
claimed in the previous section, and thus Hp is enhanced 
by a factor of \/3 compared with that in the case with 
inversion symmetry. 

4-3 Nuclear relaxation rate 1/T\ 

One of powerful experimental probes for superconduct- 
ing states is the NMR measurement of the nuclear relax- 
ation rate 1/Xi. For conventional superconductors with 
the s-wave symmetry, \jT\ exhibit a coherence peak just 
below T c , and an exponential decay at low temperatures. 
For unconventional superconductors with an inversion 
center such as High-T c cuprates, and many heavy fermion 
systems, there is no coherence peak, because of the ab- 
sence of the coherence factor and the suppressed singular- 
ity of the density of states, and also a power law decay at 
low temperatures appears. In the case without inversion 
symmetry, a quite different behavior is theoretically ex- 
pected, as pointed out by the present author and Hayashi 
etal. 22 > 23 

For simplicity, we neglect electron correlation effects, 
which may not change the following argument drastically. 
Then, the nuclear relaxation rate is given by 

1 



below T c is markedly enhanced in the case of the p-wave 
pairing dominated state. The implication of this result 
for CePtsSi will be discussed in section 5. 

4-4 Magnetoelectric effect in the normal state 

The existence of the asymmetric SO interaction 
aCo(k) ■ (T yields a nontrivial coupling between charge 
and spin degrees of freedom, giving rise to interesting 
electromagnetic properties. One example is the existence 
of magnetoelectric effects. Magnetoelectric effects in a 
normal metal were discussed by Levitov et al. many years 
ago within one-body approximation. 43 ' 44 

In the case of the Rashba SO interaction with 
Co(k) = (toy(k), —to x (k),0), the magnetization along the 
x-direction is generated by an electric field applied in the 
y-direction as 



Mr, 



(45) 

and conversely an AC-magnctic field gives rise to a charge 
current flow expressed as 

<W " (46) 



-2T 



xy 



dt 



cx 



The physical origin of these effects is easily understood 
as follows. As shown in Fig. 6 (a), when charge current 
flows along the y-axis, the Fermi surface is deformed, and 
because of the Rashba SO interaction, the distribution 
of spins changes, yielding the bulk magnetization along 
the x-axis. Conversely, an applied magnetic field in the 
lim -Im [r y^{Tr[— G(s m + u n , fc) — G(e m , fc')jc-direction changes the distribution of spins, which also 
■ ■ - - deforms the Fermi surface asymmetrically, leading to a 

charge current in the y-direction. 

In the case of heavy fermion systems, the magneto- 
electric effect coefficient T xy is significantly affected by 

/ds 1 rrAr . ,, 2 r , n2ii /, strong electron correlation. According to the analysis 
27r2Tc0sh 2 " ^ +[^a{£) \\, \ il L„„„A „„ 4-V.„ i;„„;j +1, . „.„ „Uj-„;„ „ 



-Tr[— F(e 



k.k 1 



,k)—F(e m ,k')}}\ 



iu. n —.u.+i5 



■it: 



with N n {e) and N a {e) defined by the retarded Green's 
functions as 



k t=± 

Na(e) = -J2J2^Fr R (e,k). 



(43) 



(44) 



x>ased on the Fermi liquid theory, we obtain a simple 
relation among T xy , the specific heat coefficient 7, and 
the resistivity p: 27 

akp 



Pb 



r (47) 
ev* F p E F p 

In general, for heavy fermion systems, the resistivity is 
given by p <~ p + const. 7 2 T 2 , with p a residual re- 
sistivity. Thus, at sufficiently low temperatures, T xy is 
enhanced by the factor 7. 

We now estimate the order of the magnitude of these 
effects. We assume that the Fermi velocity is v* F ~ 
10 5 cm/s, which corresponds to a mass enhancement 
of order ~ 100, and the SO splitting is sufficiently 
large, e.g., akp/Ep ~ 0.1. To consider the magnetiza- 
tion induced by an electric field, we assume that the 
charge current density is J ~ 1 A/cm 2 . Then, the in- 
duced magnetization is estimated as, M = T xy E y w 
p^iakp / Ep){J / ev* F ) ~ 1 Gauss, which is experimen- 
tally measurable. To evaluate the charge current induced 
by AC magnetic fields we assume that an AC mag- 
netic field B = B cos(ut) with _B ~ 100 Gauss, and 
lo ~ 100 kHz is applied, and that the normal resistivity 
is p ~ 10/xfi • cm. Then we obtain the charge current, 
J = -2T xy (dB/dt) w fj, B (ak F /E F )(dB/dt)/(ev F p) - 1 
mA/cm 2 . This magnitude is also experimentally accessi- 
ble. However, in this case, it is necessary to distinguish 



A^ a (e) gives the coherence factor, which enhances the co- 
herence peak of l/TfT. For unconventional pairing states 
with the angular momentum I > 2, the coherence fac- 
tor vanishes. A remarkable feature of 1/T\T for non- 
centrosymmetric superconductors appears in the case of 
the p-wave state. In contrast to the case with inversion 
symmetry in which the coherence factor vanishes for the 
p-wave state, the absence of inversion symmetry gives 
rise to nonzero coherence factor N a for the p-wave state. 
This is easily understood as follows. In the case of the 
pure p-wave state, the superconducting gap is given by 
A(fc) = AtC(k)-crio~y where C(k) has the p- wave symme- 
try. An important point is that the factor C(k) appears 
in eq. (44) only in the form |£(fc)|. (see eqs.(20) and 
(19)) Thus, when the density of states for the (+)-band 
is different from that for the (— )-band, A^ a (e) is nonzero. 
This result implies that the coherence peak of \jT{F just 
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(b) 




Fig. 6. (Color online) Changes of the Fermi surfaces and the spin 
configuration due to applied magnetic field (or net current flow), 
(a) For the Rashba interaction, current flow along the j/-axis 
gives rise to magnetization along the rc-axis. Conversely, an ap- 
plied magnetic field along the x-axis yields current flow along 
the jy-axis. (b) For the Dresselhaus interaction, the current flow 
is parallel to the applied magnetic field (or induced magnetiza- 
tion). A sectional view for k z = is shown. 



between the magnetoelectric effect and the usual eddy 
current induced by the time-dependent magnetic field. 
To avoid the eddy current, leads attached to the sample 
must be aligned exactly parallel to the magnetic field. 

A similar effect is also possible in cubic systems with- 
out inversion symmetry. In this case, the charge cur- 
rent (or applied electric fields) is parallel to magnetic 
fields (or magnetization) as shown in Fig. 6 (b), i.e., 
J = —2T(dB/dt). The coefficient T is approximately 
given by eq. (47). 

4-5 Magnetoelectric effect in the superconducting state 
Magnetoelectric effects in the superconducting state 
analogous to those in the normal state are possible. 
These phenomena were theoretically predicted by Edel- 
stein, and extensively discussed by Yip and the present 
author. 1,2 ' 6 ' 20,22 ' 27 In the absence of inversion symme- 
try, supercurrent can be induced by the Zeeman mag- 
netic field, which may be called "paramagnetic super- 
current". 1 ' 6 Conversely, bulk magnetization is raised by 
supercurrent flow. 2 The physical origin of these effects 
is also the asymmetric deformation of the Fermi surface 
due to an applied magnetic field as in the case of the 
normal state. However, in contrast to the analogous ef- 
fects in the normal state, the magnetoelectric effects in 
the superconducting state are dissipationless. As a re- 
sult, to detect the paramagnetic supercurrent flow ex- 
perimentally, one needs to prepare a system in which the 
dissipationless bulk current is possible. One example of 
such a system may be realized by attaching leads made 
of superconductors to the sample. 45 

To explain the above-mentioned effects in more detail, 
we first use the Ginzburg-Landau (GL) theory, and later, 
we will present microscopic analysis. The GL free energy 
for superconductors without inversion symmetry was de- 
rived by Edelstein, Samokhin, and Kaur et al., 3 ' 16,19 
which reads, 

1 

(HI - - 

2' 



F s - F„ 



a|*| 2 + ?l*| 4 



2m 



2en x 

B 2 

' 8tt 



X^B 2 



(48) 



where a = a (T - T c0 ), = -HV^ - 2eA^/c, is 
a vector potential, B — V x A, M is the magnetiza- 
tion density, and n s is the superfluid density. The fourth 
term of eq. (48) with the coefficient JC^ stems from the 
parity-breaking SO interaction, and is the origin of mag- 
netoelectric effects. Differentiating the free energy with 
respect to A and B, we obtain the following relations 
for the supercurrent density 3 s and the magnetization 
density M. 



Jf, = J,, 



Mf, 



(49) 
(50) 



Here, J^ la is the conventional diamagnetic supercurrent 
given by Jf a = (W M - 2eA M /c)/(2eA). <f> is the phase 



of the order parameter ^ . A 



2e 2 |*| 2 /m. Mff e is 



a magnetization due to the conventional Zeeman effect. 



Also, we have put |\&| 



The second term on the 



right-hand side of eq. (49) is the paramagnetic supercur- 
rent, and the first term on the right-hand side of eq. (50) 
is the magnetization induced by supercurrent flow. In 
the case of the Rashba interaction, K,^ is nonzero only 
for (jj,, v) = (x, y) or (y, x). In the case of cubic systems, 



JC8 n 



As was pointed out by Yip, in the case of the Rashba 
interaction, the paramagnetic supercurrent is partially 
canceled with the magnetization current J M — cV x 
M. 20 To observe this, using eqs.(49), (50), and the rela- 



tion V x J L 



-_B/cA, we write down the total current, 



r s + J M = J° 



3 T 



cV x M Zc 

dia a Tdia 



> V »y % (51) 

The last term on the right-hand side of eq. (51) is the 
paramagnetic supercurrent. In the complete Meissner 
state and in the thermodynamic limit, this term van- 
ishes, and thus there is no paramagnetic supercurrent. 
Yip pointed out that because of this cancellation, the 
penetration depth is symmetric under the transforma- 
tion z — ► — z. 20 However, in finite systems, or in the 
mixed state, the last term of eq. (51) gives nonzero con- 
tributions to the magnetoelectric effect. 

In the case of cubic systems, the cancellation be- 
tween the paramagnetic supercurrent and the magne- 
tization current is perfect. This is easily seen as fol- 
lows. Since in cubic systems /C M „ = K,8^ v , we obtain 

js = jdia + KB M _ £ AJ dia + M Zo 0j jM _ 

cV x M = -JCB + cV x M Zco , and as a result, the para- 
magnetic current K,B cancels exactly in the total current 
J s + J . There is no paramagnetic supercurrent in this 
case. 

The magnetoelectric effect coefficient JC^ can be ob- 
tained from microscopic calculations as in the case of the 
normal state. According to the Fermi liquid theory ex- 
plained in section 2, the formula of JC yx , which includes 
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electron correlation effects exactly, is given by 22, 27 
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on the basis of the GL free energy (48). From the deriva- 
tive of eq. (48) with respect to 'J, one finds that a stable 
solution is given by *(i?) = ^ a (R)e~ i( l R with 19 

fe) q>1 ~ h^ s 



l -B v . 



(54) 
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2T 
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where vo yT = d{ek+Ta\to\) / dk y . In the case with a spher- 
ical Fermi surface, up to the first order in akp/Ep, the 
magnetoelectric coefficient is simplified as 
__ e^i B n s a 



8^*25*. ' (53) 
where n s is the superfluid density. Note that JC yx is am- 
plified by the mass enhancement factor l/zk T , when the 
Wilson ratio is nearly equal to 2, as realized in typical 
heavy fermion systems. This feature is in contrast to the 
electron correlation effect on the diamagnetic supercur- 
rcnt which is suppressed by the factor Zk T - As a result, 
the magnetoelectric effect in the superconducting state 
is more strongly enhanced in heavy fermion systems with 
large effective mass than in weakly correlated metals. 

We now discuss the feasibility of experimental obser- 
vations of these effects. Assuming that akp/Ep ~ 0.1, 
the electron density n <~ 10 22 cm~ 3 , the mass en- 
hancement factor 1/z <~ 100, and v s /v* F <~ A/E F <~ 
0.01, we estimate the magnitude of the bulk magneti- 
zation induced by the supercurrent as M — /CAJ dla w 



Here ^o(R) is the solution of the GL free energy without 
the parity-breaking term (the fourth term of eq.(48)). 
This corresponds to a Fulde-Ferrel state with a spatial 
modulation of the order parameter phase. It is noted that 
(5$he period of the modulation decreases as <~ z 2 T , as 
the effective electron mass l/zk T increases. When l/q^ 
is smaller than the spacing between vortices, the helical 
vortex phase governs the upper critical field. The critical 
temperature at finite B in this phase is derived by using 
the standard method, 16,19 



fj, B n(ak F /E F )(v s /v* F )/(8n i 



0.1 Gauss. The exper- 



imental detection of this internal field may be possi- 
ble. Under the above conditions, the magnitude of the 
paramagnetic supercurrent is also accessible by conven- 
tional experimental measurements. It should be empha- 
sized again that to detect paramagnetic supercurrent, 
one needs to prepare a circuit in which bulk current flow 
without dissipation is possible. Also, in the mixed state, 
vortices should be pinned by impurities to suppress dis- 
sipation due to the flux flow, which may be induced by 
supercurrent, as demonstrated by Oka et al. 26 

4-6 Helical vortex state (Fulde-Ferrel state) 

In the previous section, we have seen that the asym- 
metric deformation of the Fermi surface due to an ap- 
plied magnetic field may give rise to the paramagnetic 
supercurrent, if the dissipationless bulk current flow is 
permitted. Conversely, when bulk current flow is forbid- 
den, as in the case of an isolated system, the asymmetric 
deformation of the Fermi surface stabilizes an inhomoge- 
neous superconducting state with a finite phase gradient. 
This possibility was first pointed out by Kaur et al. and 
Samokhin. 16, 19 They argued that in the Rashba case, 
an in-plane magnetic field applied to an isolated system 
gives rise to a Fulde-Ferrel-like state, in which Cooper 
pairs have a center of mass momentum and the phase of 
the order parameter is spatially modulated. 16, 19 Kaur et 
al. called this novel phase the helical vortex state. Fol- 
lowing their analysis, we would like to consider this phase 
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Thus, the transition temperature (or the upper critical 
field) is increased by the inversion-symmetry breaking 
term of the free energy (48). We would like to empha- 
size that the increase in T C (B) is drastically amplified 
by electron correlation effects, since the last term of the 
right-hand side of (55) is enhanced by the factor l/z| r , 
while the second term corresponding to the conventional 
orbital depairing effect is suppressed by the factor Zk T - 
Therefore, the realization of the helical vortex state is 
feasible in heavy fermion systems. To estimate the order 
of the magnitude of this effect, let us assume that, as 
in the case of CePtsSi, the ratio of the SO splitting to 
the Fermi energy ak^jEp is <~ 0.1, the mass enhance- 
ment factor 1/zkr is of order ~ 100, the Fermi velocity 
is v F <~ 10 5 cm/s, the upper critical field is H C 2 ~ 4 
T, and the applied magnetic field B is close to H c2 . 
Then, the period of the spatial modulation is estimated 
as 1/q ~ [i B B/(l6ir 3 v F z) • (akp/Ep) ~ 10~ 6 cm, which 
is compatible with the inter-vortex distance. In this situ- 
ation, the increase in T c due to the Helical vortex phase 
is almost on the same order as the decrease due to the 
orbital depairing effect, and plays a crucial role in the 
vicinity of the upper critical field. 

4-7 Anomalous Hall effect 

In general, the SO interaction can yield the anoma- 
lous Hall effect, as clarified by Karplus-Luttinger many 
years ago. 46 The anomalous Hall effect stems from the 
anomalous velocity which is caused by the SO interac- 
tion. In the case of the Rashba SO interaction, this mech- 
anism can be schematically illustrated as shown in Fig. 7. 
When an electric field applied along the y-axis induces 
current flow along this direction, the Fermi surfaces are 
deformed asymmetrically, and thus the distribution of 
spins, which is constrained by the Rashba SO interac- 
tion, is changed as shown in Fig. 7(a). Then, a magnetic 
field applied along the 2-axis gives rise to a torque which 
rotates spins around the z-axis as depicted in Fig. 7(b). 
Since the Rashba SO interaction forces the Fermi mo- 
mentum to be perpendicular to spins, the asymmetric 
Fermi surfaces are rotated on the xy-plane. As a conse- 
quence, the net current along the x-axis occurs as shown 
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in Fig. 7(c). In the case of the Rashba interaction, the 
anomalous velocity va = a\7to(k) ■ (n x er) is perpen- 
dicular to the z-axis. Thus, the anomalous Hall effect is 
possible only for magnetic fields along the z-axis. Ac- 
cording to the analysis based on the Fermi liquid theory, 
the anomalous Hall conductivity in the normal state is 
expressed as 27 



(a) 



(b) ® H z ^ 



Re a 



AHE 

xy 



e 2 MB 



EE 

r=± k 



2a|Ret( £ * T ,fe)|3 



*{dkJxTdk y toy - dkJyrdkytQx)- (56) 

This quantity is enhanced by the factor A sz , which is 
equal to the enhancement factor of the van-Vleck-like 
spin susceptibility (30). For heavy fermion systems, this 
factor is of the same order as the mass enhancement fac- 
tor l/zkr ~ 100 ~ 1000, and thus, the anomalous Hall 
effect is significantly large. For instance, let us assume 
the resistivity p ~ 10 /if2 • cm, the mass enhancement 
factor l/zkr ~ 100, the Fermi velocity v* F ~ 10 5 cm/s, 
and the carrier density n ~ 10 22 cm -3 . Then, the ratio of 
cr mf E ~ e 2 fiBB/ (h?v* F z) to the normal Hall conductivity 
cr^ HE is estimated as <r^ E /a^ y HE ~ 40. The anomalous 
Hall effect dominates over the normal Hall effect. 

An analogous effect of heat current is also possible. 
The anomalous Hall conductivity for heat current is ex- 
pressed as 27 



.AHE 



where is equal to the conductivity tensor er M „ and 



= I(£(2)_\^£(l)£(0)-l L (lh 



(57) 



- (m)AHE 
J xy 

H z 



r=± k 



2a|Re t{e* kT ,k)\ 3 



(dkjxrdkytoy - dk x t VT dk t{§} 



with m — 1,2. 

Note that in eqs. (56) and (58), electrons away from the 
Fermi surface give dominant contributions to the anoma- 
lous Hall conductivity. This feature is in accordance with 
the fact that the magnetic response against the magnetic 
field along the z-axis is governed by the van-Vleck-like 
term. This observation leads us to an interesting implica- 
tion for the superconducting state. In the superconduct- 
ing state, the Hall effect for heat current is possible at 
finite temperatures, and when the superconducting gap 
is much smaller than the size of the SO splitting, the 
thermal anomalous Hall conductivity is not affected by 
the superconducting transition. Furthermore, even in the 
limit of zero temperature, K xy ^ E /(TH Z ) is nonzero, and 
behaves as in the normal state, even though the quasi- 
particle density is vanishingly small. The experimental 
detection of this effect is an intriguing future issue. 

5. Pairing State Realized in CePtaSi 

In recent years, several heavy fermion superconduc- 
tors without inversion symmetry have been discovered. 
In particular, extensive experimental studies on the su- 
perconducting state in CePtaSi have revealed several im- 
portant aspects of this system. In this section, we would 




Ex Jx 



Fig. 7. (Color online) Schematic explanation for anomalous Hall 
effect, (a) An applied electric field along the y-axis changes the 
distribution of spin, (b) Torque produced by magnetic field along 
the z-axis rotates spins, (c) Because of the Rashba SO inter- 
action, the Fermi surface is deformed, leading to a transverse 
current along the x-axis. 



like to discuss about possible pairing states realized in 
CePtsSi, comparing the recent experimental data with 
the theoretical understanding of noncentrosymmetric su- 
perconductors. 

As mentioned in the previous sections, the d- vector of 
the spin triplet component is determined by the parity- 
breaking SO interaction. In the case of CePt 3 Si which has 
the crystal structure with the C^ v symmetry, the possi- 
ble form of the parity- breaking SO interaction aCo{k)-cr 
was clarified by Samokhin from a general group theoreti- 
cal argument. 21 According to his analysis, there are three 
constraint conditions that determine the form of Co(k): 
(i) Co(—k) = —Co(k). (ii) Co(k) satisfies the symmetry 
of the point group, (iii) Co(k) is a pseudovector which 
changes its sign under the transformation (k x , k y ,k z ) — > 
(k y ,k x ,k z ), {x,y,z) — > (y,x,z), since it is a quantum 
mechanical average of the operator VV^ x fc. 21,48 Then, 
the general form of Co(k) for C^ v symmetry is given by 
Co(k) = {ky,—k x ,cok x k y k z (k x — ky)). 21 For periodic sys- 
tems, fc M and fc 2 are, respectively, replaced by sin k^ and 



cosfc M . The coefficient Co of the z-component depends 
on microscopic detail of the system. The case of cq — 
corresponds to the Rashba interaction. In this case, the 
d- vector proportional to (k y , —k x ,0) gives the highest T c , 
as clarified by Frigeri et al. 15 In the case of Co ^ 0, the 
situation is more complicated, because the magnitude of 
the pairing interaction for the channel corresponding to 
the basis function k x k y k z {k x — ky) (the A 2 representation 
of Ci v ) is generally different from that for the p-wave 
channel corresponding to the E representation of C^ v , 
and thus d ^ A t (k)£ (k), leading to depairing effects, 
as mentioned in section 2. The pairing state in this case 
is determined by the competition between the strength 
of the pairing interaction and the pair-breaking effect 
due to the SO interaction. In the following, we assume 
that Co is sufficiently small for simplicity, and that the 
Rashba type interaction is applicable to this system, i.e., 
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One remarkable experimental observation for CePt 3 Si 
which is relevant to the pairing state was obtained from 
the NMR measurements of 1/Ti by Yogi et al. 11,12 Ac- 
cording to their data, 1/Ti exhibit a small coherence 
peak just below T Cl indicating a fully gapped state. An- 
other important experimental observation was obtained 
from the thermal transport measurements by Izawa et 
al. 13 They found that at low temperatures the thermal 
conductivity exhibits the Volovik effect, indicating the 
existence of line nodes of the superconducting gap. Also, 
the results of the penetration depth measurement, 14 and 
the specific heat measurement 47 support the existence of 
the line nodes at sufficiently low temperatures. Thus, for 
the clarification of the pairing state, it is crucial to rec- 
oncile the existence of the coherence peak of 1/Ti with 
the line node of the superconducting gap. In general, the 
coherence peak of 1/Ti is due to the coherence factor 
and the singularity of the density of states at the gap 
edge. Even for unconventional superconductors with line 
nodes, the singularity of the density of states exists, and 
may gives rise to a small coherence peak. However, in 
strongly correlated electron systems, because of the large 
quasiparticle damping, the singularity is smeared out, 
and as a result, no coherence peak is observed for all un- 
conventional superconductors in heavy fermion systems 
except CePtaSi. In CcPt 3 Si, the quasiparticle damping 
is quite large as indicated by the resistivity data in the 
normal state. Thus, the existence of a coherence peak 
cannot be explained by the singularity of the density 
of states. Also, it is very unlikely that the s-wave pair- 
ing is realized in such heavy fermion systems with the 
large onsite repulsion. Therefore, we need to give some 
specific reason for this remarkable feature. As argued in 
section 4.3, in noncentrosymmetric superconductors, the 
coherence peak of 1/Ti is enhanced even for the p-wave 
pairing dominated states, because of the non-vanishing 
coherence factor. Thus, from the argument presented in 
section 3, one possible candidate for the pairing state re- 
alized in CePtaSi is the s + p wave state, in which the 
p-wave pairing is predominant. The present author and 
Hayashi et al., independently proposed this pairing state. 
Furthermore, Hayashi et al. pointed out that in the s+p 
state, when the magnitude of the s-wave gap A s is close 
to that of the p-wave gap A p , a line node of the excitation 
gap appears, because A_ = A s — A p sin# can be zero for 
a specific 9, where 9 is the azimuthal angle. Hayashi et al. 
attributed the line nodes observed at low temperatures in 
CePtaSi to this mechanism. 23 This scenario is tempting, 
since all the experimental observations are explained only 
by tuning one parameter A s /A p . However, as claimed 
above, the substantial weight of the s-wave admixture is 
a bit unlikely in heavy fermion systems in which there is 
a strong onsite repulsion. The present author proposed 
another scenario in which the coupling with an antifcr- 
romagnetic (AF) order that coexists with the supercon- 
ductivity in CePtaSi gives rise to nodal structures of 
the excitation gap. 49 The drastic change in the excita- 
tion spectrum due to the coexistence of the AF order 
is not inherent in noncentrosymmetric superconductors, 
but was pointed out for centrosymmetric conventional 



superconductors previously. 50 In the case of CePtaSi, the 
AF order is characterized by the ordering Q- vector along 
the (001) direction and the staggered moment parallel 
to the a6-plane. 9 We demonstrated that for sufficiently 
large magnitudes of the staggered moment, a line node 
structure appears in the vicinity of the magnetic Bril- 
louin zone even when the superconducting gap itself has 
no nodes. Assuming the p-w&vc pairing state coexisting 
with the AF order, we calculated the energy dependence 
of the density of states D{e). The results are shown in 
Fig. 8. Here, tuq is the magnitude of the staggered mo- 
ment. In the calculation, we assumed the Rashba type 
SO interaction, and the d-vector is d — (k y ,—k x ,0). At 
sufficiently low energies, D(e) cx e holds; a characteris- 
tic of line nodes. In this calculation, it is assumed that 
the staggered moment is aligned parallel to the (100) 
direction. In Fig. 9, we show a schematic figure of the 
line node structure of the Fermi surface. As a matter of 
fact, this node structure is not a true node, but a mini- 
mum of the excitation gap. Thus, the coherence factor is 
not suppressed by this node structure, which is in accor- 
dance with the existence of the coherence peak of 1/Ti. 
It is noted that the node structure has Ci v symmetry in 
momentum space because the direction of the staggered 
moment breaks C\ v symmetry. The experimental detec- 
tion of this two-fold symmetry is a crucial test for this 
scenario. 

The possibility of the realization of the s+p wave 
state should be examined by microscopic calculations. 
For CePt 3 Si that has a tetragonal crystal structure with 
Ci v symmetry, basis functions for the s-wave states cor- 
responding to the irreducible representations for A\ g are 
1, cosfca; + cosfc y , and cosfcz, as clarified by Sergienko 
and Curnoe. 18 Thus, a possible candidate for the gap 
function is given by eq.(14) with A s (fc) = A^ + 
Ai 1} (cosfc^ + cosfc y ) + Ai 2) cosfc z , A t (k) = A| 0) + 
A^(cos fc^ + cos k y ) + A[ 2 ^ cosfcz, and the Rashba type 
C{k) = (sinfc y , — sinfc x , 0). Even for this gap function, 
since the isotropic s-wave component A^ is generally 
nonzero, and gives rise to the coupling with the strong 
repulsive onsite interaction, its stability should be exam- 
ined by model calculations. Unfortunately, it is difficult 
to construct a reliable microscopic model for CePtaSi at 
the present moment, since experimental studies of mi- 
croscopic electronic structure such as the de Haas-van 
Alphcn effect have not yet succeeded in detecting the 
Fermi surface of heavy electrons in CePtaSi. 35 Very re- 
cently, Yanase and Sigrist have carried out microscopic 
calculations using the Fermi surface obtained by the LDA 
calculation, and found that the admixture of a p-wave 
state and an extended s-wave state is stabilized. 51 

6. Concluding Remarks 

In this review, we survey our theoretical understand- 
ing of the unique properties of noncentrosymmetric su- 
perconductors. One of the main conclusions is that the 
remarkable effects due to parity-violation are strongly 
enhanced by electron correlation, and thus may play cru- 
cial roles in heavy fermion systems. From this point of 
view, the recently discovered superconductors without 
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Fig. 8. (Color online) Density of states for p-wave pairing state 
plotted against excitation energy e. The upper panel is on a linear 
scale. The lower panel is a log-log plot. 




Fig. 9. (Color online) Schematic figure of line node structure on 
the Fermi surface. The depth of the blue color in the vicinity of 
the magnetic Brillouin zone k z ± 7r/2 indicates the depth of the 
node of the excitation gap. 



inversion symmetry such as CePtaSi, CeRhSi3, CeIrSi3, 
and Ulr are quite fascinating systems. It is highly desir- 
able to reveal the interesting features inherent in noncen- 
trosymmetric superconductors experimentally in these 
novel systems in the near future. 

We should note that this review, by no means, en- 
compasses all topics on noncentrosymmetric supercon- 
ductors. In particular, we could not present arguments 
on issues on the mixed state in detail. Semiclassical ap- 
proaches to these subjects have been explored by Hayashi 
et al. and Nagai et al. 24,52 In the case of A/ Ego *C 1 as 
realized in all heavy fermion superconductors without in- 
version symmetry, the pairing states on different Fermi 
surfaces can be approximately separated even under ap- 
plied magnetic fields, and thus the theoretical treatment 
is much simplified. However, when A/ Ego ~ 1 5 with fi- 
nite magnetic fields, the pairing between the different 
Fermi surfaces can not be negligible, and presumably 
yields drastic changes of low-energy properties. As a mat- 
ter of fact, as pointed out recently by Eremin and Annctt, 
the Zeeman field destroys the line node structure due 
to the admixture of the s-wave pairing and the p-wave 
pairing considered by Hayashi et al. 23 and mentioned in 



section 5, leading to a full-gap state. 53 The theoretical 
clarification of orbital effects in such cases is an interest- 
ing future issue. 
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